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Abstract. The chemodynamical evolution of spherical multi-component self-gravitating models for isolated dwarf 
galaxies is studied. We compare their evolution with and without feedback effects from star formation processes. 
We find that initially cuspy dark matter profiles flatten with time as a result of star formation, without any 
special tuning conditions. Thus the seemingly flattened profiles found in many dwarfs do not contradict the cuspy 
profiles predicted by cosmological models. We also calculate the chemical evolution of stars and gas, to permit 
comparisons with observational data. 
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1. Introduction 

Dwarf galaxie s are the most common galaxies in the 
1 Universe (e.g. iMarzke fc da Costal (|1997l )) and the Local 
Group sample is an excellent labo ratory where t o study 
these systems and their properties ( Mateo ( 1998 ): Grebell 
ljl999h l. Our interest in these structures is related to many 
important mutually related topics: the growth of struc- 
ture, the nature of the dark matter and dynamical in- 
teractions. In the last decades, new observational data 
have become available for many dwarf galaxies (in par- 
ticular in our Local Group) especially regarding the ve- 
locity dispersion profiles of stars and the rotation curves 
of gas, providing clean measures of the dynamical mass 
at all radii where baryonic matter can be found. This is 
gradually permitting a deeper investigation on the equi- 
librium and stability of the dark matter profiles. Such 
studies revealed a discrepancy between the prediction of 
the KCDM models by numerical simulations, which re- 
quire power law den sity profiles w i th cusp s p oc r~ 7 where 
1 or 1 .5 (e.g. iMoore et al. ( 1999a ): Navarro et al, 



(|2002l ): lGentile et alJ (|2005h : lde Blokl (120051) 1 . Other prob- 
lems strictly related to the nature of the dark matter and 
its clustering properties are, e.g. , the missing s atellite s 
problem (e.; 



Moore et all (Il999al): iKlvpin et al.1 (|l999al) 



(119971 120041) 1. and th e more flatten e d profiles derived 
from observations (e.g. iBurkert ( 1995 ): de Blok fc Bosmal 



Moore et all (|1999bl ): IKlvpin et al.1 (|1999b|) 1. the trivial- 
ity of the halos for cl usters of galaxies in ferred from gravi- 
tational lensing (e.g. iTvson et al. I(ll998l)1 and the angul ar 
momentum problem (e.g. INavarro fc Steinmetj (|2000l) 1. 
Despite the many successes of the KCDM paradigm, 
a large number of open questions remain. Among the 
large body of literature on the nature of dark matter 
and tes ts of differe n t formulations of gravity, th e work 
of, e.g.. iBode et all (l200ll) : ISpergel fc Steinhardtl (l2000h 
and iBekensteinl (|2004f ) provide classical results. Our goal 
in this paper is to extend these studies by examining the 
interplay between dynamics and star formation in the evo- 
lution of dwarf galaxies. Although including star forma- 
tion makes the modeling more complex, it also allows for 
comparisons with dwarf galaxies in the Local Group. 

Chemodynamical modeling of this type has been pro- 
ductive in the past, particularly in investigations of the 
effects of winds from hi gh-mass stars an d supe r nova explo- 
sions m the ISM (e.g. IPelupessv et~all (|2004l) : ISlvz et al 
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(|2005l) : Ide Avillez fc Breitschwerdtl (|2005l )). Given the 
small dynamical mass inferred for dwarf galaxies, it 
is surprising that the ISM can remain bound long 
enough to permit star formati on episodes that l ast fo r 
a few Gyr (as fo und by , e.g., Babusiaux et al.l (200a); 
Carnev fc Seiteen (Il986l) : iDolphinl <|2002h : iMateol Jl998), 



1 



Grebel et al.l (|2003l )) giving a typical baryonic matte 
binding energy of ~ 10 53 erg with a total expected number 
of SN events of 10 3 4- 10 4 . 

We will focus our attention in this paper on the dwarf 
spheroidal (dSph) galaxies. These systems are the least 
massive galaxies known but their velocity dispersions 
suggest a mass-to-light ratio of up to 1OOM /L , 
making them some of the m ost dark matte r dom inated 
objects in the universe (e.g., Gilmore et al.l (2007)). For 
these galaxies, the star formation history and chemical 
composition are very environment dependent, making 
the realization of a consistent theory of their origin 
and evolution that combines dynamics, star formation 
and chemical enrichment more difficult. High-resolution 
spectroscopy of s everal dSphs sh o ws a large spread in 



metallicity (e.g. 
by photometr ic studies 



Aparicio et al.l (|200ll ); iBellazzini et al 



Shetrone et al.l (l200lh) corr ob orate d 



Harbeck et al.l 
(I2002K 



(l200ll) : 

Several 

authors suggested that the ISM of d Sph systems could b e 
entirely remov e d by S N e explosions (iDekel fc Silkl (11986); 



Hensler et all (120041) : iMori et al.1 (120021 120041. Il997l): 



Murakami fc Babul ( 1999 ) ; Lanfranchi fe Matteuccil 



( 2004: iMac Low fc FerraTaN l 999)1 



Mashchcnko et al.1 (1200a) discuss how star formation 



processes can successfully act as flattener for the central 
density cusps, within very short timescales (more details 
will be given in the following). iRead fe Gilmord (j2005l) 
obtain the sam e result as effect of an external impulsive 
mass loss event. Gnedin fc Zhao (12002 ) investigate the in- 
fluence of winds. Navarro et al.l ( 19961 ) were one of the first 
to point out the possibility that feedback mechanisms can 
turn the central dark-matter cusp into a cored one. 

We start our analysis by demonstrating the gravita- 
tional stability of the isolated model. Then, once the stel- 
lar processes are included, we combine the dynamical and 
chemical analysis. In Section[2]we outline our research ap- 
proach. In Section[3]we set up the initial model, in Section 
[3] we present the evolution of the models with and with- 
out stellar processes, in Section [5] we add the chemistry 
analysis and in Section [51 we discuss our results. 

2. The modeling approach 

In this paper we are trying to address the problem of the 
dark matter profiles in the dwarf galaxies introduced in the 
previous Section. Despite the complicated processes act- 
ing on the true physical systems, different tools have been 



1 Hereafter we identify baryonic matter with the stellar com- 
ponent; the dark matter will generally be called not-baryonic 
matter even if its nature could be baryonic as in the case 
of baryonic-dark-matter objects such as MACHOs or brown 
dwarfs. 



developed and presented in the literature to follow the 
chemodynamical evolution of an isolated spherical multi- 
component self-gravitating model of a dwarf galaxy sys- 
tem. Briefly, we will base our arguments on a chemody- 
namical code. Our code can follow the time evolution of 
dynamically interacting smoothed-particles that represent 
stars, gas and dark matter components of a dwarf galaxy. 
The code takes stellar evolution and chemical enrichment 
of the stellar populations and of the gas into acc ount. A 
full de s cription of the code u sed can be found in iBerczik 
(1999); ISpurzem et al. (2009), to which we refer the reader 
for details on codification and the implementation. We will 
use observational constraints in order to tune our models 
and simulations. We will proceed as follows: 

1. We first realize multi-component self- gravitating mod- 
els of a dwarf galaxy composed of a cuspy dark matter 
density profile, a stellar density distribution and a gas 
component. 

2. We evolve these models in isolation to check their grav- 
itational stability. 

3. We activate the star formation (SF) processes to sec 
the changes in the dark matter profiles (if any) and 
vary the physical parameters. 

4. We compare the results to understand the interplay of 
SF processes and dynamics. 

5. We turn off the stellar evolution processes to check 
again the dynamical equilibrium of the resulting sys- 
tem. 

Our goal is to explore how star formation influences 
the properties of the dark matter density profiles by chang- 
ing the baryonic density profiles. 

The same approach will be tested on an initially cored 
dark matted density profile (see later sections.) 

Our focus on isolated systems is the first hypothesis 
that makes our modeling not immediately comparable 
with observations since completely isolated galaxies prob- 
ably do not exist in nature. 

Our first hypothesis leads to the following corollaries: 

1. We can assume a '7-profile' (see below) that represents 
quite well both cored and cusped dark matter profiles 
in the inner part of the dwarf galaxies. The inner slope 
of a 7-profile has the advantage of depending on just 
one single parameter, i.e. 7. For the gas and stellar 
density profiles adopted in the starting models we can 
again assume gamma profiles as good representation 
of the projected luminosity profile of the dwar f galax - 
ies in the de-Vaucoul eurs form (e.g. Dehnen ( 19931) : 
Tremaine et~aTl (|l994T )) 



The outer part of the true dwarf galaxy is not expected 
to drop off as the 7-models. However, we are consid- 
ering a synthetic spherical isolated model, so that for 
Newton's first theorem at any given radius f the dy- 
namics within f are not affected by the mass distri- 
bution exterior to f. As a consequence we expect a 
priori no differences in the evolution of isolated mod- 
els regardless of the slope by which our dark matter 
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profile falls off, in particular in the outskirts of the 
galaxy whenever spherical symmetry is assumed. This 
permits us to work consistently with the 7-models. 

We point out that these two corollaries are not true 
if the first hypothesis is not assumed. Nonetheless in the 
literature often different dark matter profiles are used, 
especially from the family of the King profiles and the 
Plummer spheres. But, even if the gravitational force at a 
given radius is not dependent on the matter distribution 
outside of that radius, this is not true for the tidal forces. 
Therefore, we expect that in real systems, the tidal influ- 
ence of external galaxies is relevant for the determination 
of both the shape of the DM density profile as well as its 
inner slope (except in the unlikely case of a dwarf galaxy 
evolved at rest exactly at the center of a spherical clus- 
ter of galaxies) . The tidal effects are expected to act both 
by truncating the dark matter profile by tidal compres- 
sion and by stretching the outer part of the dwarfs when 
evolved as satellite systems of a primary. Thus, changes 
in the dark matter profiles in the directions of the lead- 
ing and trailing tails are also expected, or shallower pro- 
files in the outer part of the dwarfs have to be used if we 
want to introduc e our isolated system in the cosmological 
framework, e.g., iNavarro et al. I (|l997t) (hereafter NFW). 
However, that is not our objective here (but see the com- 
ment in section 2]). 



3. Isolated galaxy 

3.1. The model 

Our goal is to provide a self-gravitating three-component 
model of a dwarf galaxy. The evolution proceeds initially 
without any astrophysical process other than gravity. This 
is done in order to check the consistency of the model and 
to permit us to isolate the different effects in the evolution 
of a true astrophysical system due to gravitation or stellar 
evolution. 

In order to choose realistic parameters for our model 
of a dwarf galaxy, we look at a system that in the liter- 
ature has been subject of extensive investigation, i.e. the 
Carina dwarf galaxy. Our target is not to realize a spe- 
cific model of the Carina dwarf galaxy, but first to use 
observationally motivated input physics and then extend 
our analysis to the dwar f galaxies class of systems (but 
see 



Piatek et al] (|2004h ) . Dynamical considerati ons are avail- 



Pasc tto et al. The wealth of literature data 



available on this galaxy, starts already in 1983 with studies 
of the stellar populations of Carina ( Mould fc Aaronsonl 
(1983)). Color magnitu de diagrams were used to derive the 
star formation histo ry ( Hurlev-Keller et al.l|l998t iMighelll 



19971 I1992L fl 990 



story 

M 



Rizzi et alJ 



Morielli et all l2006bllal . I2004bllal ld. 
2004 120031) 1 Spectros copic mea- 



2003allbT 

surements y ielded chemical abunda nces, e.g..lKoch et al 



(|2006l 120081 ). the M/L ratio, e.g., iGilmore et all (120071) 



and t he 
( 2003h 



dark matter 



Gilmore et all 



distri bution, e.g., 
(l2007h. Proper motion 



Havashi et al 



nations were obtained by iPiatek et all (|200 



determi- 
ne also 



able f or the orbit from timing argum ents in|Munoz et al 



(|20081 ) and for Carina's tidal tails in iMunoz et all (|2006f ) 
which ma ke this system well-sui ted for further orbit inves- 
tigations ([Pasetto et al.l (|2009ah ) 

We are adopting the following density distribution pro- 
file: 

I \ 3-7 Mr s 



4tt r i (j + 



From the analysi s of the CMD of Carina (e.g. 
Hurlev-Keller et al. (1998)) we infer that the oldest stel- 



lar population contributes no more t han 15% of the stel - 
lar mass. Moreover from the work of Mateo et al. I (Il993h 
we get a mass to light ratio in the V band of (77) y ~ 

23 (j^j . Thus from Table 3 of lHavashi et all (|2003h 



we get L = 0.43 • 1O 6 L0, which leads to a stellar mass 
estimate of M « 23 • 0.43 • 1O 6 M = 9.89 ■ 1O 6 M . If 
we now assume that the mass lost by tid al interactions 



could reach an order of 95%, following e.g.. lHavashi et al 



(|2003l ) we get M Car = 2 ■ 10 8 M Q . From this initial mass 



only 1/50 is assumed to be observable. This baryonic ma- 
terial consists of 15% stars, 85% gas. We thus adopt a 
starting value of Af gas = ^^^^ = 3.3 • 10 6 A/ Q and 
Af star = ^^-^ = 5.9 • 10 5 M Q for our reference model 
These numbers can be c onsidered as starting values for 
an orbiting dwarf galaxy ( Pasetto et al 2009al) . but here 
we will work only on an isolated system, thus for us the 
gas amount retained by the dwarf galaxy after its oldest 
stellar population is formed will be considered as a free 
parameter. This approach is assumed because the present- 
day dSphs contain little or no gas, possibly because a large 
fraction of gas was expelled by the initial SNII component, 
by tidal interaction with a host galaxy or in general pho- 
toevaporation, rampressure, by star formation processes. 
Depending on its own individual evolution and on its 
orbital parameters, each dwarf ultimately evolves under 
different conditions depending on its orbital peri-center 
passages and a general study cannot be included in the 
present work. 

The initial temperature for the gas can be inferred 
assuming a spherical collapse model for the initial dark 
matter model and a matter power spectrum P(k) at red- 
shift zero compatible with studies from the SDSS, e.g., 
iTegmark et al. ( 2004) or work on the analysi s of the 
Lyman-a forest, e.g., Gnedin fc Hamilton (|2002f) . In this 
case the typical RMS internal velocity of a halo within 
M < 1O 9 M is less than 20km-s _1 . This implies that e.g., 
hydrogen with a sound speed higher than 20km ■s~ 1 has a 
velocity c s greater than the escape velocity of the poten- 
tial well in which it has to collapse. This gas is hence not 
in the condition to be retained by the dwarf. Therefore we 
will assume an initial sound speed in hydrogen of roughly 
around c s = 10km -s -1 corresponding to a temperature of 
T ^ 10000°K. 

The evolution of the reference model presented here 
will be used as normalization case for the realization of a 
Carina-like dwarf galaxy model that is discussed in more 
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detail in Paper II (jPasetto et al.l (|2009a)). Carina is one of 
the few dSph galaxies in the Local Group that shows long- 
lasting star formation and pronounced intermediate-age 
populations. This mean Carina resembles to much more 
isolated, low-mass dwarf irregulars and is one of the rea- 
sons why we chose it as a reference for the approximate 
masses of the different galactic components. We empha- 
size that it is not the intent of this paper to reproduce 
the star formation history of the actual orbiting Carina 
dSph. That is the subject of a separate paper. In fact, in 
our simulations we explore a larger range of stellar, gas 
and dark matter masses (Table [!} in order to bracket th e 
classica l dwar f spheroidal galaxies given in iMateol ([19981 ) ; 
Grebell (Il997h . 

Another parameter on which Eqn. ((T|) depends, is the 
scale radiu s r c . It is inte r esting to point out that in the 
analysis of lHavashi et al. (2003) of the possible extension 
of the dark matter halo, they deduce that the tidal ra- 
dius r t of the assumed NFW shaped halos can greatly 
exceed the surface brightness cut-off i? co , by an order of 
r t > 20R CO and that the rate of mass loss can vary from 
5% to 95% depending on the specific orbit chosen. We 
proceed by tuning the outer profiles for our dwarf model 
as follows. We can express the radius containing 90% of 
the mass, rgo, as function of the scale radius for the pro- 
files of Eqn. ([T]) with 7 = 3/2 by solving the simple Eqn. 
4tt J r9 ° p (r) r 2 dr = 9 / 1Q M, which yelds: 



r c = — [27 + 10 • 30 1 / 3 + 3 • 30 2 / 3 ] »•.,„ 



Total Mass [M©] 


from 


to 


7 




A/atar 


5.1x10" 


1.0x10 s 


7 = 


= 1.5 


Afdark 


7.1xl0 7 


1.9xl0 9 


7 e [0.5,2] 


M gas 


5.5xl0 4 


6.9xl0 6 


7 = 


= 1.5 



For o ur model we adopt r t = rgo (jlrwin &: Hatzidimitriou 
19951) and we assume an initial tidal radius 20 times 
as large as the present one (jHavashi et al.l 120031 ). With 
these values the stellar and gas profiles have rgo = 
(0.5 x 20) kpc = 10 kpc, leading to a core radius of 
r C) T — 0.72 kpc at the starting time To of our simula- 
tion^. 

Since we are not interested in the external radial de- 
pendence of the profile, as previously remarked, we could 
produce the same outer radial dependence for the dark 
matter profile by modifying the scale length in order to 
obtain the same r go , radius containing the 90% of the mass 
as used for the baryonic profile. However the approxima- 
tion adopted for massive particles with Plummer spheres 
imposes a natural limit on the inner resolution radius we 
are going to analyze. The assumed smoothing length for 
our model (see below) is roughly sum — 0.120pc depend- 
ing on the specific model. We will not attempt an analysis 
within the limiting radius r; = 2edm i n the following. 
Nonetheless we need to be able to discriminate between a 
dark matter profile with 7 = 2 or 7 = for radii r > rj . 7- 
models with an assumed total mass of Mum — 2 • 10 9 A/ Q 
are almost indiscernible for r G [n,^9o] for rgo = 20 kpc 
at our particle resolution. 



Table 1. Range of masses within which the simulations 
are run. The scale lengths refer to r 90 as discussed in the 
text. The values listed here are our starting values at the 
initial time T = Tq 



4. Evolution 

We will focus our attention here only on the properties of 
interest for dwarf galaxy systems, for which the following 
results hold. The role of the triaxiality and the extension 
to the more general class of elliptical galaxies is not in- 
vestigated in the present paper. A table with the range of 
masses considered here is presented in Table Q] Within 
this range we mostly sampled models with an jt = 
[Afdark/Af s tar]T e [16,688]. This range of ratios partially 
excludes very massive 'dwarfs', e.g., Large-Magellanic- 
Cloud-sized proto-sy stems (unrealistic for t he scale radius 
adopted here but see lPasetto et al. I (|2009bl) ) and we com- 



The definition of To is provided in the following section. 



pletely excluded system with initial M s tar+Af gas > -Afdark- 
The age of the Universe is 13.7 Gyr, or if expressed 
as lookback time, T, n , = —13.7 Gyr. We start our sim- 
ulations much later, namely at a time when the dwarf 
galaxies we want to analyze have already formed their 
old populations. Our final goal is to study the influence of 
star formation on the dark matter density profiles. In oder 
to isolate more clearly this effect we need to neglect the 
Hubble flow, and all the heating/cooling processes that 
influence the star formation at hight redshift. Hence, by 
starting at lookback time of To = — 9Gyr we can follow 
the phase-space description of our models in a simpler 
way and implement the star formation processes and their 
feedback on the interstellar medium with higher numer- 
ical resolution. Nevertheless, this approach forces us to 
introduce a preexisting stellar population in order to take 
into account the gas evolution from Tj n j to To. In this ini- 
tial time the AC DM paradigm predicts that the dwarf 
galaxies are the first structures to form as scale invariant 
perturbations having masses greater than the Jeans mass 
of isothermal perturbations at the time of recombination 
(masses of = lO 5 "^ 7 Mq). Hence, in this scenario we need 
to consider the time to violently relax the dark halo com- 
ponent within which the proto-cloud will congregate and 
produce the first stars. It is clear that these complicated 
phases, treated in cosmological simulations, make it more 
difficult to clearly evidence and disentangle the single stel- 
lar evolution influence on the whole gravitational potential 
of a single system because two fundamental conditions we 
want to use are missing: equilibrium and isolation. After 
this initial phase and after an equilibrium configuration 
has been reached, say at To = — 9Gyr, we expect that an 
old stellar population has already formed. 
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4.1. The evolution without the star formation 
processes. 

Once the model has been realized we evolve the sys- 
tem in isolation in order to test the equilibrium of the 
system's density profiles. The critical parameter is of 
course the smoothing length of the particles, although in 
the treecodes the gravitational force estimation is largely 
performed through cell-cell interactions. Different stud- 
ies in the literature can be found on th e best-fit tunin g 
of the smoothing length p a ramete r (e.g. Defined (|200lh : 
lAthanassoula et ah (2000); Zhan ( 2006 ) and references 
therein). We mention here that different codes require to 
be tuned on different softening lengths to grant the sta- 
bility of the starting model depending on the implemen- 
tation of the code. The initial conditions evolved on a 
special purposed hardware (GRAPE) reach stability with 
e S 0.12 kpc. 

As can be seen in Figure [1] (left panel), stabil- 
ity is achieved within the resolution used here. The 
figure represents a simulation with 35000 particles 
(n gas ,n sta r,n-dark) = (1.0,0.5,2.0) x 10 4 . The same sta- 
bility is obtained with slightly different values for the 
smoothing length e of the dark matter, stars, and gas 
with a larger number of particles (n gas , n star , "-dark) = 
(2.0, 1.0, 4.0) x 10 4 and (n gas , n star , n dark ) = (1-0, 0.5, 2.0) x 
10 5 for a total amount of 350000 particles in the beginning 
of the highest resolution simulation. Even if the higher res- 
olution would permit us to analyze regions closest to the 
core of the dwarf galaxy, we will see that our chosen dis- 
tance of 0.5 kpc from the core to which we limit our con- 
sideration due to the smoothing length resolution of the 
dark matter particles, is enough to confirm our results. In 
Figure [T] (left panel) , the overlapping lines represent the 
evolution of the system at T = To, i.e. 9 Gyr ago, and 
after AT = 9 Gyr, i.e. now (T = T + 9 = -9 + 9 = 
Gyr) , respectively. The density profile of the model hardly 
changes with time as is desired, both for dark matter pro- 
file as for the stellar component, thus confirming that the 
initial conditions of our model are reasonable. 



4.2. The evolution with the star formation processes. 

Once the equilibrium of the density profile has been con- 
firmed under the effect of gravity alone, we activate the 
star formation processes. 

The understanding of the star formation processes was 
a primary task of the Astro nomy during a l l the secon d 



half of the last century (e.g., iLarsonl (|l969lh Isilkl ^98^)) 



but, despite the great improvement of knowledge acquired, 
the micro-physics description still represents a challenge 
of the modern astrophysics. Great improvement has been 
achieved within the framework of the Smooth Particles 



Hydrodynamics s t ar for mation algorithms (jKatzl (|1992l ): 
Navarro fc White! (|1993l )) that we adopt in our descriptio n 
in its more updated recipe as explained in lBerczik (1999). 

We emphasized again that we start our simulation af- 
ter the old stellar populations have already formed. For 



instance, if it took the dwarf galaxy w 3 or 4 Gyr to ac- 
crete gas, to have it collapse and to experience its first 
measurable episode of star formation, then we assume this 
equilibrium stage as our initial time of evolution starting 
from which we analyze the system (T = To), 3 to 4 Gyrs 
after the Big Bang. 

The assumed dynamical equilibrium of the model is 
as previously described in Section 14.11 for the old stellar 
component, the gas and the dark matter, under the effect 
of gravity alone. The system can now evolve from an initial 
time T) until the present time. Since we are interested in 
systems with cuspy dark matter density profiles, we only 
present their evolution here. Further considerations on flat 
dark matter profile systems can be found in Section [5] 

Depending on the initial amount of gas and on the total 
mass of the system, the general behavior of the system is 
to expel part of the gas due to the effect of SNII that grad- 
ually eject the gas to the outer part of the galaxy. Some 
of the gas is nevertheless conserved in the inner part of 
the system thanks to the steep potential well. The dwarf 
galaxy starts as an extremely dark matter dominated sys- 
tem especially in the inner part of the galaxy where the 
density profile is very cuspy (M/L > 200 within 5 kpc 
from the center) . In this sense, the steepness of the dark 
matter profile seems to be a necessity in the initial phase 
of the dwarf galaxy where just a few percent of the stellar 
population are formed after the primordial collapse. There 
arc at least three reasons for this: 



1. 



If the initial dark matter density profile is not cuspy 
but flat, the possibility for the dwarf galaxy to survive 
several pericent er passages as expected for a dwarf like 
Sagittarius (e.g. lHelmi fc White ( 200ll) ) becomes more 
problematic. 

If the initial dark matter profile is not cuspy but 
flat, the efficiency of the S NII to expel the gas is ex - 



pected to be higher (e .g. iBradamante et al" (jl99 



Read fc Gilmore (120051^ thus increasing the amount 



of primordial gas necessary to reproduce tidally trig- 
gered bursts of star formation in galaxies as may be 
the case in, e.g., Carina (see Paper II). However, this 
result is strictly dependent on the initial environment 
in which the gas collapses and we start our simulations 
only after a time To where all these effects have already 
taken place. 

Cosmolog ical N-body simulatio ns predict cuspy pro- 



files (e.g. iNavarro et al 



(11997^ with an inner slope 
proportional to r « r . The asymptotic behavior for 
the 7 models in the inner zones predicts, in the ac- 
ceptable range studied here, r oc r -7 , including, by 
the way, the NFW inner slope. 

It is difficult to constrain this residual fraction of gas. 
It seems to be a free parameter that can play a role in 
th e periodic bursts of st ar formation recently suggested 
bv lPelupessv et all (|2004l) or may play a role in the dwarf 
evolving in an orbit where the SFR can be activated by 
strong tidal interactions. 
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Fig. 1. (left panel) Density profile temporal evolution for the model evolved in isolation under the action of gravity 
alone. The blue lines represent the dark matter profile at the beginning (solid) of the simulation and at the end 
(dashed). The red lines show the stellar profile. No apparent evolution of the dark density profile is seen confirming 
the equilibrium of the model simulated. The baryonic component remains exactly the same as the starting profile, 
hence the red dashed line is not visible here since it coincides with the solid red line, (right panel) Density temporal 
evolution for the same model of the right panel evolved in isolation and considering the most relevant star formation 
processes. The blue solid line represents the dark matter profile at the beginning of the simulation T = Tq. The final 
profile is the blue dashed line. The same line types are used for the stellar profile (red lines). Intermediate thinner 
lines refer to intermediate evolutionary stages: T = T + 2 .9 Gyr and T = T h + 6.1 Gyr The graph scales are chosen 
to be directly comparable with similar results obtained in Mashchenko et al. ( 20061 ). 



As we can see from Figure [T] for the reference model 
of Section [3J we experience a stellar density profile evo- 
lution with the same time scale, t s , as the dark matter 
density profile time scale = t s = 9 Gyr. The system 
evolves slowly over timescales of a few Gyr modifying its 
stellar profile thanks to stellar evolution: new stars are 
born and die interacting with the small fraction of gas to 
reach each time new configurations of dynamical equilib- 
nu n§ We point out here how not all the literature recipes 
for the SF are able to follow efficiently the gravitational 
changes of the baryonic component. The code really has 
to 'make' the stars such that the number of particles in- 
creases regardless of the substantially increasing compu- 
tational time required. Onl y these prescriptions, s ee e.g . 
Figure H and iBerczikl dl9 99l): iBerczik fc Kravchukl (|2000h; 
Berczik fc PetrovT ( 200ll ). together with the correct grav- 
itational treatment (here provided by a devoted GRAPE 
hardware) permit one to follow correctly this slow change 
of the gravitational potential. In Figure[2]we plot the range 



3 This is obtained whenever the cell opening parameter of the 
treecode keeps a conservative value, say 6 « 0.6° in order to 
grant a better energy conservation, or when a more extensive 
use of the G RAPE hardware is adopted (e.g., IBerczik et al.l 
(2006, 2007)) 



of particles used to test our results and their time evolu- 
tion. The thick lines refer to the reference models and the 
thinner ones to the test models as defined in Section 14.11 
The test model were developed in order to confirm the 
particle number independence of our results. The slight 
variation on the number of dark matter particles is due 
to a small evaporation in the outer part of the dwarf that 
occurs at distances out of 10 kpc from the center of the 
system. The number of particles for the stellar and gas 
components changes depending on the star formation pro- 
cesses (Section 15. ip as well as on the evaporation in the 
outer part of all the components (gas, stars and also dark 
matter). The similarity of the evolution any trend is re- 
markably good from the reference model up to the test 
model with the highest particle number. 



Once new stars are added to the stellar component, 
they reshape the mass distribution (the density profile) 
and as a consequence the dark matter reacts with the 
same scale time t s = passing thought several states 
of dynamical equilibrium. This can be proven by simply 
stopping all the star formation processes and evolving the 
new equilibrium configuration under the influence of grav- 
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Fig. 2. Variation of the number of particles for the refer- 
ence models and higher resolution test models. Black solid 
lines refer to the total number of particles, green dashed 
lines refer to the gas, red dotted lines refer the stars and 
blue dot-dashed lines refer to the dark matter. The refer- 
ence model is marked with a thicker line. 

ity alone. The phase-space configuration seems to be sta- 
ble with respect to the simple dynamical evolution^. 

We observed that the inner flattening of the dark mat- 
ter halo is accompanied by a simultaneous decrease of the 
density profile of the baryonic component (Fig [1} . That is 
simply an indication that the final density profile has to 
depend at least on two parameters (other than the total 
mass or the central density). As an example, the family of 
7 models or King models permit this kind of metamorpho- 
sis: we can obtain this result by simply fixing the radius 

containing 90% of the mass rgo = r s (-^) 3-7 — 1 and 
expressing in this way the scale radius as a function of rgo 
and 7. At fixed rgo and total mass M in Eqn. ([1}, one can 
easily obtain central densities that differ by a few orders 
of magnitude. 

More intuitive is the explanation of the evolution of 
the stellar density profile because it changes both central 
density and scale length due to the increasing number of 
particles and mass, while nonetheless roughly preserving 
the projected surface brightness. 

In Figure[3]we plotted the rotation curve v c — J r^| of 
the dark matter component at the starting time T = Tq 
and at the end of the simulation T = Tq + 9 Gyr (the 
intermediate-time snapshots have been omitted for the 
sake of clarity). The general trend expected as a conse- 



4 This might seem a numerical test on the stability of the 
equilibrium configurations but it is not because it is model de- 
pendent. The treecode approach prefers the cell-cell force com- 
putation for distant interactions by reducing the role played by 
the softening length (which is, on the other side, dependent on 
the growing number of particles used in the systems) . The sta- 
bility of the equilibrium configurations will not be subject of 
further investigation in this paper. 
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Fig. 3. Rotation curve for the dark matter profile (and the 
stellar component) of the reference model. The solid line 
represents the starting profile at T = T and the dashed 
line the profile after T = Tq + 9 Gyrs. The same colors 
and line types are used as in the previous figures. 

quence of the flattening of the profile is exactly the de- 
crease of the central part of the rotation curve, as shown 
in Fig. [3] We can however observe the small overlap of 
the rotation curve, which in a two parameter density- 
potential pair is an indication that the matter distribution 
at T = Tq + 9 Gyr is slightly more concentrated. A fit of 
the dark matter profile shows that the outer homoeoids 
for the dark matter profile have slightly contracted. The 
trend is generally present in our simulations even though 
the effect here is minimal. Nevertheless, this is a further 
indication of the fact that a good dark matter density 
profile has to be at least a 2-parameter model. 

The cumulative mass density profile M (r) = 
4.7T J Q r r' 2 p(r') dr' shows this trend as well (Fig. |4|. The 
theoretical expectation for the curves in Fig. 2] indicates 
that the flatter density profile curve asymptotically con- 
verges under the cuspier density curve profile (once we 
assume the same total mass). Considering the asymptotic 
behavior, the overlapping part of the final cumulative mass 
profile is an indication of a slightly reduced final scale ra- 
dius: the same amount of matter is reached in inner ho- 
moeoids and starting from a central smaller value as shown 
in Fig. HI 

Before discussing the kinematics we point out the lim- 
its of our simulations. Our finding that a two parameter 
dark matter density profile is favored, is based on 

— the flattening of the density profile at a constant total 
amount of mass. 

— the analysis of the rotation curve and cumulative mass 
profile. 

But only the first point is stringent evidence for the ne- 
cessity of the two free parameters model in the dark mat- 
ter component (in addition to the total mass), while the 
second point comes with an intrinsic theoretical bias and 
can not be used to claim the necessity to add a further 
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Fig. 4. (upper panel) The cumulative mass profile for the 
dark component (blue) and the stellar component (red). 
The starting profile is at T = Tq (solid lines) and mostly 
lower than the profile for time T = To + 9 Gyr (dashed 
lines), (lower panel) A zoom on the dark matter compo- 
nent shows how the lines of the initial and final stage of 
the evolution cross. Sec text for details. 



parameter to the NFW profile. The low reactivity of the 
outer homoeoids to the changes of the inner mass distri- 
bution is expected as direct consequence of Newton's sec- 
ond theorem. This theorem describes how the total force 
acting on the matter distribution out to, e.g., rgo, is the 
same as it would be if the shells' matter were concentrated 
into a point at the center of the dwarf galaxy, despite 
the different radial distribution it shows in spherical ap- 
proximation. Thus the finding that the stable dark matter 
profile has to be at least a two-parameter density profile 
is just an indication that comes from the baryonic/non- 
baryonic interaction analysis and does not necessarily rule 
out th e one-parameter density profiles as in lNavarro et all 
( 19971) . Moreover as already explained in Section a 
model evolved in isolation is expected to show a sharper 
decline of the external profile as compared to the one re- 
sulting from the cosmological collapse of a cube (as we ef- 
fectively see by comparing a 7- model with a NFW model). 

To investigate these processes further in detail we 
present in the following also some chemical considerations 
that will help us to improve the present theory on dwarf 
galaxy evolution. The idea to include some chemical con- 
siderations, emerges from the analysis of the anisotropy 
parameter, defined as f3 = 1 — s^f, where ay is the radial 



velocity dispersion and of = c^+cr^ is the tangential com- 
ponent in a suitable spherical coordinate reference system 
(O, r, tp). It is evident from this definition, and from as- 
suming that, as for every reasonable Galactic potential 
of > of, that the anisotropy parameter (3 6 [— 00, 1], 
where /3 — — 00, corresponds to a system dominated by 
pure circular orbits ay = 0, in contrast to a system with 
(3 = 1 with pure radial orbits. As seen in Fig. [5] the sys- 
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Fig. 5. Anisotropy parameter (light green) and radial ve- 
locity dispersion (red) for the stellar component. The left 
scale on the y axis is for the anisotropy parameter, the 
right scale for the radial velocity dispersion parameter. 
The thick solid lines show the starting stellar profile, here 
at To = 196 Myr, and the dotted lines show the resulting 
profile at T = T Q + 9Gyr. 



tern develops, within the timescale of the flattening of 
the dark matter cusp, a zone dominated by circular or- 
bits, roughly out to a radius ru m — 4 kpc in radial dis- 
tance. The kinematics in the outer part of the galaxy are 
left unchanged in this isolated modeH and the radial dis- 
persion velocity profile remains basically unchanged be- 
yond r > ru m . Similarly the anisotropy parameter ap- 
proaches unity at larger radii, getting closer to unity after 
T > Tu m = 4kpc. This indicates how circular orbits be- 
come more and more relevant when the dark matter cusp 
becomes larger, offering a plateau of uniform density and 
potential, where we expect an easier mixing of the stellar 
component. This will leave some observable traces in the 
radial gradient of the chemical composition that we now 
are going to explain. 



5. Chemo-dynamical evolution 

5.1. Star formation recipes 

The coding of the star formation processes is not the 
subject of our analysis. A more ext e nded d escription can 
be foun d in th e works by Berczik J 19991); Berczik et al 



(l2003al lbL l200l : lBerczik k Petrovl (l2o6l¥lBerczikll200d) . 



Nevertheless the robustness of our evolution has been 
tested against the different star formation criteria that 
are discussed in the literature. We found that there are 
several physical conditions that we have to take into ac- 



5 Any further investigation of the link between the scale 
length of the starting model r c and this limiting distance is 
not the subject of this paper. Nevertheless when evolving this 
isolated system in an orbit around a host system, we expect 
that the external tidal force can perturb any relation that may 
exist between this limiting radius ru m and the scale length r c . 
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count in order to obtain the evolution of the dark matter 
profiles: 

— A fluid element has to be instable, e.g., it has to sat- 
isfy the Jeans instability criteria T soun d > Tff between 
sound speed timescale and free fall timescale. 

— For the gas particles t coo i < r//, the cooling timescale 
has to be lower than the free fall timescale (this is also 
often expressed in terms of the over-density criterion 
Pi > Pcriti i- e -j the density of the gas particle i, p%, has 
to be higher than a critical density p cr it)- 

— The efficiency of the energy that heats the ISM in a 
dwarf galaxy can be reduced to a few percent e = 0.03, 
say erjsNiji where tjsn is the standard energy injec- 
tion per SN explosion. This is a result that was derived 
by the chemical modelin g approach already presen ted 



, x 10" 
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>• 0.6 



to 



0.4 



0.2 




in the literature, by, e.g.. iBradamante et alJ (|1998f l 



— A radius limit of the sphere of influence of a SNII 
event can be assumed when the sound speed of the 
front shock is cooled at the sound speed of the ISM. 
Simulations bv lCioffi fc Shulll (|l99lh have shown that 
a reasonable value is of the order of a few hundred 
parsec. 

These are sufficient although not necessary conditions. 
There are other studies that correctly reproduce the 
observed dynamical effects, e.g., when using dynamical 
times, Tdy n — y/3n/16Gp instead of the fr e e fall time 
in the first criterion (e.g., iNavarro fc White! (|l993h h or 
when using convergent fluxes with the crit erion of the di- 
vergen ce of the velocity V • v < (e 



Carraro et al 



T Q+ T [Gyr] 

Fig. 6. Evolution of the star formation rate with time for 
the reference model. The total amount of mass (the sum 
under the area of the histograms), can change by an or- 
der of magnitude depending on the free parameter of the 
amount of gas that the dwarf galaxy can retain, e.g., in 
our case the initial amount of gas. Nevertheless it is always 
less than a few percent of solar masses per year. The total 
absence of spurious initial bursts of star formation is also 
a confirmation of the good level of dynamical equilibrium 
obtained in Section EOT) 



led to gas ejection to the outer regions of the dwar f 



Recchi et alJ (|200ll )). The first two criteria are uni- 



(e-g- 

versally accepted and do not require any particular fine 
tuningin our dwarf galaxy simulations to reproduce our 
resultqj (see Fig. The last two criteria are not tuning 
parameters. They are results that were already presented 
in the literature many years ago. They basically act as SN- 
'shock absorbers' in the dwarf galaxy environment. Thus 
the result we have presented for the dark matter evolu- 
tion is not dependent on any physical tuned mechanism 
acting on the baryonic component, nor does it require any 
'special' mechanism to explain the change of the baryonic 
density profile and in turn the resulting gravitational ef- 
fect it produces (the conversion of a cuspy dark matter 
profile into a flat DM profile). 

5.2. Chemical Evolution 

In our simulations we assume that the first episode 
of star formation and its associated SN explosions 



galaxy as described in, e .g.. iMac Low fc Ferrara 



(|1998)). or even applymg higher values of SN efficiency |L an f ranch i fc Matteuccil (|2004h : iMarcolini et al 



1999); 



2006) 



The only slight modification of the Jeans instability cri- 
terion is in the case of the presence of magnetic fields. In this 
situation also the first criterion can be slightly modified in fa- 
vor of a critical mass depending on the Alfven and the sound 
speeds. Despite their relevance in the collapse and fragmen- 
tation of cold clouds, magnetic considerations have not been 
investigated in our simulations. 



This is not a direct result of our simulation because, as 
mentioned, we consider the evolution of a dwarf galaxy 
after its initial star formation episode which gives rise its 
old stellar population allowing us to better control the gas 
present at the initial stage T — Tq of our simulations. 

In Fig. [7] we consider the metallicity evolution with 
time, without considering its spatial dependence (see Fig. 

by plotting the trend of iron enrichment. Newly 
formed stars peak above the mean value of iron, ( [Fe / H] ) 
which as consequence increases with time and is strongly 
dependent on the initial chemical composition assumed 
for the ISM and for the stellar component. We are actu- 
ally not interested in reproducing any specific observed 
trend, but we want to understand the expected rate of en- 
richment. Starting with a rather metal-poor population 
({[Fe/H]) = —3 dex), we are able to see the general 
rate of enrichment independent from the specific initial 
conditions that differ from galaxy to galaxy in the true 
Universe. It is evident that a mean range of ([Fe/H]) E 
[—1.5, —0.5] can be covered in 9 Gyr of isolated evolution 
provided that the system has an initial range of rather low 
metallicity and then undergoes rapid early enrichment in 
the initial star formation episodes th at need to have taken 
place d uring the first, say, 3 Gyr (e.g.lMac Low fc Ferrara 
(|l999h : lLanfranchi fc Matteuccil (|2004fi ;i. We note that ob- 
servationally, such a rapid initial increase does seem to 
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T 0+ T [Gyr] 

Fig. 7. The evolution of the stellar [Fe/H] over time. The 
contour plot shows the metallicity spread distribution with 
time. At the starting time To the dwarf starts with a metal 
poor distribution associated with an old stellar population 
formed prior to the beginning of our simulation. The color 
scale is normalized to the highest number of stars (that is 
not constant in our simulations!) and to the peak of the 
metallicity per bin. In this way the plot better demostrates 
the increase of the metallicity with the time and its dis- 
tribution per temporal bin (due to the gradual filling of 
the distribution tails by SNIa events). 



take place, see, e.g., iKoch et al. but we need 

to remark that a direct comparison with any real dwarf 
galaxy can not be carried out here. With the trend here 
evidenced we want to reproduce an 'upper-limit' to the 
observed chemical enrichment. Thus the results here pre- 
sented are reference model that can be used to impose an 
upper limit constraint. 

In Fig. [5] we plot the number of SN events and PN as 
a function of time. This helps to explain how the spread 
in the previous figures originate s. In line with the argu - 



ments of different authors (e.g. iMarcolini et al.l ((2006)) 



we see how the SNIa, mainly responsible for the iron en- 
richment of the ISM, operate to fill gradually the tail of 
the (([Fe/H]) — t) diagram, even though their number is 
only a few percent of the SNII. This can indeed be de- 
duced from the Figure [5] where the number evolution of 
SN events is presented as a function of time. Depending 
on the chemical enrichment we assumed for the oldest 
stellar population, at the end of the 8 to 9 Gyr of evo- 
lution we can obtain a mean present day metallicity of 
(\Fe/H\) = —0.5 integra t ed along the radial direction (c fr . 



Lanfranchi fe Mattcucci 



IMarcolini et al.1 (|2006)) 



The initial dominance of SNII is partially suppressed 
in this model thanks to the ad-hoc condition avoids 
the initial burst of star formation (see Fig and the 
SNIa rate is decrea s ing w ith time as, e.g., predicted in 
Matteucci fc Recchi (2001)). The reason is that in the col- 
lapse of the proto-cloud the star formation is enhanced by 
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Fig. 8. Cumulative number of SNII events in red, SNIa 
in magenta and planetary nebulae in gray. The number 
of SNII is almost constant due to the assumption of con- 
tinuous star formation after the initial collapse which is 
assumed to to have occurred prior to time To, by creating 
the old stellar population of our dwarf galaxy. 
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Fig. 9. Here the SNIa rate is presented. We find agreement 



with the literature prediction (e.g. iGreggio fc Renzini 
1 1983 )) of a decreasing S NIa rate with time (see also 
Matteucci fc Recchil (j200ll) ) 



the increasing density (SFR oc p 1 ^ 2 ) thus favoring the SN 
events, while in our simulation this phase is neglected and 
a starting Population II sample is always initially present 
at T = T . 

Prior to the onset of contributions of SNIa, stellar 
metallicities ([Fe/H]) tend to increase while the [a/Fe] 
ratio remains rougly constant. This leads to a fiat, plateau- 
like feat ure at low metallicities when pl ot ting [Fe/H] vs 
a/Fe ] (jLanfranchi fc Matteuccil <|2004h : IMarcolini et al 



<|200fih ). Since our simulations only start after the forma- 
tion of the old stellar population, our Fig.[TU]docs not show 
such a plateau. Instead we see the characteristic decline 
of [a/Fe] with increasing [Fe/H] due to the combined ac- 
tion of Fe production in SNIa and a element production 
in SNII. 
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Fig. 10. The evolution of the a elements (represented by 
oxygen) and iron for different initial gas fractions. We 
show the evolution of the mean [Fe/H] and [a/Fe] in 
the gaseous and stellar components. The characteristic flat 
trend of [a/Fe] at low metallicities is not seen here since 
our simulations start after the first significant episode of 
star formation. The total mass of the gas is incremented 
from dark to lighter red: Af gas = 5.81 • 10 5 M o , M„ 



7.31 • 10 5 M Q , M, 



'0, 

3.23 • 10 e M, 







6.55 ■ 1O 5 M , M g 
and M gas = 6.46 • 10 6 M Q . The supersolar evolution comes 



from a test model with M. 



12.8 • 10* M m inside the 



usual scale length of r Stgas — 0.51 kpc (see Section^ 



Moreover, differences in the initial gas fractions affect 
the metallicities reached during the first episodes of star 
formation. Lower gas fractions imply less enrichment. Fig. 
[TOl shows results for different initial gas fractions of 15%, 
30% and 45% as an example. These result in differences 
in the amount of gas that the galaxy is able to retain at a 
fixed dark matter amount. Of course, the true behaviour of 
a dwarf galaxy depends on its intrinsic properties but also 
on the tidal interaction that may activate SF processes. 



5.3. Radial chemical gradients 

We now turn to the radial dependence of the chemical 
composition. The spatial distribution of stellar popula- 
tions of different ages and metallicities is a matter of 
debate in different studies (see below). Here we want 
to show what is expected from the simple evolution of 
our isolated synthetic model. Observationally one tends 
to find gradients such that the younger, and/or more 
met al-rich populations ar e more centrally c oncentrated 



(e.g. lHarbeck et al.l ([20011 ): iKoch et all (120061) 1. The gen 



eral expectation is that the SNe generate local enrichment 
at the beginning of the evolution, and that then diffusion 
tends to erase any radial gradient with time. This affects 
both iron and oxygen but on different time scales. While 
oxygen is produced by SNII and so distributed across the 
entire dwarf galaxy quite fast, iron is expected to be more 




Fig. 11. Radial dependence of the mean [Fe/H] value. 
[Fe/H] is averaged along spherical shells. In red we have 
plotted the stellar metallicity radial gradient from T = To 
(thick line) to T = T + 9Gyr (the dashed line) with steps 
AT = 1 Gyr . 



concentrated in the centr al regions where the S NIa are 
more densely located (e.g. Marcolini et al. (2006)). As re- 
sult the gradient is expected to be stronger for iron and 
longer lasting, while smaller and shorter in duration for 
oxygen. 

In this analysis we suggest a possible link between the 
timescale of the change of the dark matter density pro- 
files and of the formation of chemical gradients which also 
involve the evolution of the stellar orbit families. If we con- 
sider Fig. [TT] together with the evolution of the anisotropy 
parameter, in Fig.[SJ we see how the progressive growth of 
the number of stars and the resulting changes of the stellar 
density profile cause a change in the overall gravitational 
potential, leading to a new equilibrium configuration for 
the dark matter, which changes its profile from cuspy to 
flat. For the stellar component we can study the evolution 
of the orbits, which helps us in interpreting the gradient 
and the link with the dark matter. As can be seen in Fig. 
1111 in the first few Gyr, the gradient is present only in the 
central zones, where the star formation is more intense due 
to the higher density. The spatial extent of the gradient 
encompasses up to 1 kpc after 1 Gyr, up to 2 kpc after 3 
Gyr etc. This effect is initially due to the higher fraction of 
radial orbits (Fig [5]) which mix the inner and outer regions 
of the system with high efficiency As time passes the new 
density profile generated by the newly born stars gives rise 
to a growing fraction of circular orbits (preferred when a 
system is slowly evolving, because they are closer to the 
equilibrium and at lower energy) . This can also be inferred 
from the slow decrease of the radial velocity dispersion 
which leads to a decrease of the anisotropy parameter. 
This is happening within the limiting radius ru m , within 
which the chemical gradient and the anisotropy parameter 
develop for the 9 Gyr of our simulation. Once the fraction 
of circular orbits is increased (corresponding to the de- 
crease of the anisotropy parameter seen in Fig. [5]), then 
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other hand wc can not reach the resolution of the chemical 
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Fig. 12. Shells of isovelocity for the sound speed of the 
gas component. The sketch shows a cut tilted by 45 de- 
grees w.r.t. the vertical axis Z in order to evidence the 
overlapping shells of isovelocity (see text for details) . 



the mixing efficiency decreases, because there are fewer 
orbits that extend to the galaxy outskirts. This, together 
with the gradually decrease of the gas fraction, acts to 
slow the radial expansion of the zones where the chemical 
gradient can be evidenced: as seen in Fig. [11] the lines for 
T=7 Gyr, T=8 Gyr, T=9 Gyr lie closer and closer to the 
dotted line (T = T + 9 Gyr) for r < 5 kpc. Evidence of 
such a radial gradient is already present in the literature 
on N-body simulations, see e.g., [Ch iosi fc Carrarol §002) 
where the massive spectrum of the elliptical galaxy and 
the role of the stellar winds are further investigated. 

In our simulation we did not include galactic winds, 
thus excluding them as possible explanation of the chemi- 
cal mixing and consequently a faster disappearance of the 
[Fe/H] gradient. We point out that recently the same re- 
sult has been found wit h a different mechanism , but still 
avoiding stellar winds bv lMarcolini et al (2008). We agree 
with the general result presented in their paper and extend 
it by showing that the reshaping of the dark matter profile 
seems to work on the same time scale of the metallicity 
gradient formation. The overlap of the shells of the SNII 
is also evidenced in Fig. [12] In this figure an especially in- 
tense episode in the time evolution around Tq + T — 0.730 
Myr is plotted where two SN events are clearly visible (in 
red) and the remnants of at least other two other ones are 
still evident (in orange/yellow). Different spots of higher 
temperature/sound velocity are shown in the velocity sur- 
face cut arbitrary tilted by 45 degrees w.r.t. the vertical 
axis Z. This is a convenient representation of the velocity 
(or temp erature) map t hat is partially comparable with 
Fig 3 of ICioffi fc Shui] (Il99lh but here in a full 3D de- 
scription. We clearly see how two green shells of irregular 
shape (sound speed contour of 25km • s _1 ) from two dif- 
ferent SN remnants are overlapping. 

On the one hand, we are able to follow the behavior of 
the orbits and the evolution of the density profiles. On the 



codes used by iMarcolini et al] (|2006h thus preventing us 
from confirming the general efficiency of the mechanism 
invoked by these authors. Nevertheless we sustain their 
complementary picture because we also find the possibility 
of overlapping polluted regions of the SNII (see, e.g., Fig 
[T2")) . In this figure the isocontours of the sound speed veloc- 
ity of the gas component are shown. A limiting radius for 
the sphere of influence of the SNII event feedback on the 
ISM is assumed when the sound speed of the front shock 
is cooled at the so und speed of the ISM , here 11km • s 



(see Section 131 and ICioffi fc Shull (|l99lh ). This permits it 
to have overlapping phenomen a of the shock fronts in the 
dwarf galaxies as proposed by IMarcolini et al. I (l2008h . 



6. Conclusions, discussion and limits of the 
approach 

In this paper, we focused our attention on the problem of 
the dark matter density profiles by investigating the ap- 
parent incongruence between dynamical theoretical pre- 
dictions of a cuspy dark matter profile and the observa- 
tions that indicate a flatter profile. With our simple de- 
scription based on the 'classi cal' approach of smoothe d 
particles hydrodynamics (e.g. Hernauist fc Katz ( 19891 ) ) 
we have obtained the results that we briefly summarize 
here. 

The classical description of the star formation pro- 
cesses seems to be sufficient to change a cuspy profile into 
a flatter one. This is shown with the classical smoothed 
particles approach to the hydrodynamic equations of mo- 
tion completed with the star formation criteria based on 
the Jeans theory. The evolution of dark matter profiles 
could offer a natural explanation for t he appar e nt fla t 
density profiles found in observations ( Burkert (Il995l); 
Ide Blok fc Bosmal (|2002h : iGentile et all (|2005l k Ide Blok 
(l2005h once applied in a realistic context. If we intro- 
duce our theory in the context of the cosmological hierar- 
chical merging picture, where the dwarf galaxies are the 
first stellar population systems that form, the theory de- 
veloped here predicts that today the star formation pro- 
cesses should have flatted all the initially cuspy dark mat- 
ter density profiles (for the range in masses explored). Any 
dwarf galaxy that contains a stellar population older than, 
say 8 Gyr, is expected to have a flat dark matter profile, 
because either the tidal gravitational effect has changed 
the dark matter profile, or because internal star forma- 
tion processes acted on the dwarf in isolation as presented 
here. The combination of both effects is in fact expected 
to accelerate the profile flattening by tidally activating 
and enhancing the star formation processes. We point out 
here that the dark matter density profile change happens 
as soon as the star formation processes act to modify the 
gravitational equilibrium of the stellar component and, 
in turn, the gravitational equilibrium of the dark matter. 
Once the star formation efficiency is reduced, i.e. the stel- 
lar over gas fraction exceeds the smaller value predicted in 
TableHJ then the effect is minor. An initially very gas defi- 
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cient dwarf galaxy can not experience this profile change. 
On the other side, galaxies that reach dynamical equilib- 
rium with a high amount of gas still bound to the system 
(say M gas = 10 7 M Q at T > T ) do not have observational 
counterparts because they produce a dwarf system with 
super-solar ([Fe/H]) (see, e.g., FiglTU)). Nevertheless dwarf 
galaxies with an initial old stellar population of less than 
1O 4 M probably belong to a different dynamical regime 
that is not of interest for us here. 

By analyzing the time evolution of the dark matter 
profiles, we extensively investigated their dynamical prop- 
erties and evidenced the necessity for at least a two free 
parameter density model (in addition to the total mass) in 
order to correctly represent the shape resulting from the 
baryonic/non baryonic interaction. This is a result par- 
tially dependent on the isolation of the model and does 
necessary rule out NFW profiles, which are obtained only 
in a baryon-free simulations. 

In our simulations, we used an extremely strongly 
peaked cuspy DM profile in order to study the effect ex- 
pected in a case of a very deep potential well for the dark 
matter profile. As consequence these phenomena act even 
on smaller time scales for a central shallower dark matter 
density peak, or for small M/L ratios. The same effect, 
in a different cosm o logica l context, has been described by 



Mashchenko et al 



(l2006ll . although they derived a very 
short timescale for the effect (rg < 800 Myr). The most 
plausible explanation for their extremely short timescale 
Td is an 'overheating' by SN of the interstellar medium. 
The role of the SN feedback (especially the SNII sphere 
of influence and the tjsmii) are more carefully taken in 
consideration in our present study as compared to earlier 
work by other authors. This leads us in a natural way to 
a longer Td, which presents attractive advantages. In the 
context of dwarf galaxies tidally perturbed by the grav- 
itational field of a primary host (i.e., a major galaxy, or 
a cluster of galaxies) the existence of deeply embedded 
systems (e.g. Sagittarius in our MW) could become more 
problematic whenever the orbital pericenter lies close to 
the center of the hosting sys t em if r<j is as short as pre- 
dicted bv iMashchenko et al. (2006). Our models do not 
suffer this problem because the cuspy profile in our case is 
lasting several Gyr and is gradually flattened. As a result 
we are easily able to save the bu lk of the orbiting dw arf 
galaxy around a primary system (|Pasetto et al.ll2009al ). A 
second differ ence is that we treat comp letely isolated mod- 
els while the IMashchenko et ah (2006) have to work with 
interacting objects. This offers us the advantage to elimi- 
nate spurious perturbations and to have a clean determi- 
nation of the role of star formation alone. The peculiar 
orbital path of every different dwarf galaxy when evolved 
in interaction has a determining role on its SFH and on 
the density profile evolution. Our modeling of isolated sys- 
tem offers easier analysis and theoretical interpretation of 
the results. 

We tried to link the dark matter profile change with an 
observable property. A suggested indication in the paper is 
the chemical radial gradient of metallicity. We performed 



an experiment that shows how a [Fe/H] radial gradient 
may originate, in absence of galactic winds, as simply a 
result of the baryonic reshaping of the density profiles on a 
timescale compatible with the dark matter profile reshap- 
ing. Thus, in this model, the presence of a radial gradient 
is always expected in the evolved dwarf as consequence 
of a smaller anisotropy parameter and a flat dark matter 
profile shape. Finally, we want to point out a few general 
caveats: 



The fraction of gas has generally to be assumed as a 
free parameter (within the range indicated in Table 
[T|) both in isolated and orbiting dwarf galaxies. We 
expect that different contents of gas left in the pri- 
mordial dwarf galaxy can cause different star forma- 
tion histories especially when the syste m evolves in 



orbit a r ound a more massive host (e.g. iMaver et al 



( 2001bl ): IPasetto et all (|2003l l1. leading to differences 
in the baryonic and dark matter phase space re-filling 
We start with an already present stell ar population 



formed in isolation. We will show in Pasetto et al 



(2009a) ' bat the gas consumption is mostly led by the 
orbit evolution. 

Another type of observed SFH is the one presented by 
local systems as Dra and UMi. These systems show 
long-lasting early star formation episodes that led to 
considerable chemical enrichment. Obviously a sub- 
stantial amount of the gas must thus have been re- 
tained in this early period (of, say, 13-10 Gyr). When 
this type of system lies within the well of an host- 
ing system, the major role in the density profiles 
evolution is surely played by the external environ- 
ment. Even in absence of gas, the tidal stripping is 
able to smooth syste m irregularit i es and to convert 



them in dSph (e.g., IMaver et ail (|2001a|) 1. This ef- 



fect is shown to be efficient also in the MW dwarf 



companions (e.g., IPasetto et al.l (|2003h ) by morpho- 
logically evolving their density profiles. By extending 
these arguments to the more crowded cluster of galax- 
ies (as Coma or Virgo) , we enter the galaxy harassment 
regime ( Moore et al. 19961 ) that is a pure gravitational 
effect on which the literature has still put too few con- 
straints from the star formation and chemical point of 
view. 

While the trend that leads to the chemical radial 
gradient is universally evident across the range of 
masses explored (and also generally present in the 
mo re extended sample of t he elliptical galaxies, see. 
e.g. Chiosi & Garrard (1200211 its gradient slope is not. 
^E^JJu shows a strong dependence on the initial frac- 
tion of gas and stars that is not studied here. 
The possible role of the triaxiality of the system is not 
investigated here. This is expected to play a minor role 
for the range of masses here investigated. Nevertheless 
the link between the ru m , as defined in this paper, and 
the scale length of the model r c for the models here 
adopted has still to be investigated even for the more 
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general class of ellipti cal galaxies and is the subject of 
a forthcoming study ( Pasetto et al. ( 2009bl) ). 
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